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ABSTRACT 

We analyse V and //-band surface photometry of a sample of 18 Sb-Sd galaxies. Combin- 
ing high resolution HST images with ground-based NIR observations, we extract photometric 
profiles, which cover the whole disk and provide the highest possible resolution. This is the 
first photometric study of late-type spirals for which the stellar kinematics have been mea- 
sured. For 10 out of the 18 galaxies, HST data in both F160W (H) and F606W (V) are 
available, and, for those, we present colour maps and radial colour profiles at the resolution 
of the Hubble Space Telescope. 

Colours vary significantly from galaxy to galaxy, but tend to be highly homogeneous 
within each galaxy, with smooth and flat colour profiles. Some of the colour maps show jumps 
in the inner regions, likely due to dust. We determine extinction-maps in an almost model- 
independent way using the V ~ H colour map and the SAURON Mg b absorption line map 
of Ganda et al. (2007). The maps show that Ay ranges from to 2 mag, in the center from 
to 1.5 mag, in agreement with the models of Tuffs et al. (2004). 

We describe the surface brightness profiles as the superposition of an exponential disk 
and a Sersic bulge. The bulges are small (0.1-2.5 kpc), and show a shape parameter n ranging 
from w 0.7 to 3, with a mean value smaller than two: well below the value for the classical 
de Vaucouleurs bulges. Most galaxies (16 out of 18) show a central light excess above the 
Sersic fit to the bulge, which can be interpreted as a nuclear cluster, as shown by previous 
studies. We provide zero-order estimates for the magnitude of these components. We discuss 
the correlations among the structural galaxy parameters and with other relevant quantities 
such as Hubble type and stellar velocity dispersion. We compare these results with a recent 
paper by Graham & Worley (2008), who present a summary of most of the near-IR surface 
photometry of spirals in the literature. For both early and late type spirals, bulge luminosity 
strongly correlates with central velocity dispersion; at constant velocity dispersion, later-type 
bulges are larger and less dense, and have lower Sersic n values. 

Key words: galaxies: bulges - galaxies: evolution - galaxies: formation - galaxies: photometry 
- galaxies: spiral - galaxies: structure 



1 INTRODUCTION 

It has been known for many years that the inner regions of spiral 
galaxies are fundamentally different from their outer parts. The in- 
ner parts are generally rounder, redder, and show a higher fraction 
of random vs. ordered motion. Images of spirals show that less star 
formation is present in the central regions. Radial surface bright- 
ness profiles can be fitted well by a large exponential disk and a 
central, steeper component (e.g., Kent 1984). When looking at the 
stellar kinematics, one finds that v/cr in the inner parts is much 
lower than in the outer parts, where the rotation velocity is high 



and the stellar dispersion low (e.g., Noordermeer 2006). As a re- 
sult, people traditionally viewed spiral galaxies as a combination of 
a flattened disk and a spheroidal central component called 'bulge', 
which are assumed to be physically and dynamically different: in 
this picture, the disk component is rotationally supported against 
gravity, and the bulge is a hotter system, similar to an elliptical 
galaxy. To understand how spiral galaxies are built, one generally 
does a bulge-disk decomposition and studies the parameters of both 
components. If one wants to go into more details, one can go one 
step further, and also study possible components s uch as bars, spi- 
ral arms, rings, a nd inner disks (see for example Ide Jong|[l996al 
IPrieto et al.ll200ll and Erwin & Sparke 2002). 
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When doing the bulge-disk decomposition, it is imperative 
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that the decomposition is well-defined. The bulge-disk decompo- 
sition can be based e.g. on the surface brighteness profile (photo- 
metric decomposition), the stellar kinematics (kinematic decompo- 
sition) or on the two-dimensional axis ratio (morphological defini- 
tion) (see Peletier 2008). The physical interpretation of the result- 
ing parameters then depends strongly on the definition. Up to now 
most B/D decompositions have used the photometric decomposi- 
tion, using a one-dimensional azimuthally averaged surface bright- 
ness profile. Such profiles can generally be fitted well by an ex- 
ponential disk and a central Sersic r^'^" distribution (Andredakis, 
Peletier & Balcells 1995, Graham 2001). The Sersic index n has 
been found to correlate with the morphological type of the galaxy, 
and generally ranges between 1 and 4. For a comprehensive review 
about the Sersic index see Graham & Driver (2005). In recent years 
it has been discovered that the bulges with surface brightness pro- 
files that are close to exponential (i.e. Sersic index close to 1) are 
different from bulges with larger Sersic indices, which resemble 
more elliptical galaxies. Bulges with exponential profiles appear to 
contain more dust, show more recent star formation, are more flat- 
tened and more rotationally supported (Kormendy 1992, Kormendy 
& Kennicutt 2004, Fisher & Drory 2008). These objects are called 
pseudo-bulges, but one might also call them central disks. For more 
detailed examples and references we refer the reader to the compre- 
hensive review by Kormendy & Kennicutt (2004), which presents a 
complete summary of observational evidence for this kind of bulge 
and their formation via secular evolution: bulges could result from 
the evolution of disk dynamical instabilities. Numerical simula- 
tions seem in fact to suggest that the dissolution of bars inside the 
disks may trigger the formation of thr ee-dimensional stellar struc- 
tures with roughly exponential profiles jPfenniger & NormanI 199(1 
ICombes etal.lll990l iRaha et all 1 199 1. Norman, Sellwood, Hasan 
1996). It is also known that these processes are more effective in 
late-type galaxies rather than in early-types, despite the fact that 
there are examples of pseudobulges in Sa and even SO galaxies 
terwin et al 2003, K ormendv & Kennicuttll2004l) . In this paper, we 
will also adopt the photometric bulge-disk decomposition. 

Our purpose in this paper is to investigate the nature and in- 
terconnection of disks and bulges in a class of rather poorly stud- 
ied objects: late-type spiral galaxies, full of dust and star forming 
regions, and characterised by relatively low-surface brightness; in 
the last years, galaxies towards the end of the Hubble sequence 
have been the targe t for a number of HST-based imaging surveys 
dCaroUo et al1ll99llCarollo & Stiavelli 1998, Carollo, StiaveUi & 



Mack 199 8, 'Carollo| ll999i. iCarollo et al.ll200d iBoker et al.ll2002l 
iLaine et al. 2002) which revealed the presence of a variety of struc- 
tures in their inner regions: bulges, nuclear star clusters, stellar 
disks, small bars, double bars, star forming rings, whose formation 
and evolutionary patterns are not properly understood yet. Because 
of the extinction, a photometric analysis based on near-infrared data 
gives more reliable results {Ah is roughly a factor 8 lower than 
Av'- Rieke & Lebofsky 1985). Another advantage of working in the 
NIR is that the NIR light is much less affected by recent episodes 
of star formation than optical light, and traces the old stellar pop- 
ulations hosting most of the luminous mass of galaxies. For the 
analysis of HST-NICMOS data the problem is, however, that the 
field of view is very small (about 20" x 20" for NIC2, so that an 
accurate determination of the sky background is very difficult. Here 
we use data from the 2MASS survey to solve this problem. 

The sample that we study here is special, in that accurate stel- 
lar kinematic maps, stellar absorption line maps of a few strong 
lines, as well as emission line strengths of [OIII] and H/3 have been 
determined, making it possible for us to study the relation between 



structural parameters and stellar populations, the behaviour of ion- 
ized gas and the stellar potential. These maps come from SAURON 
(Bacon et al. 2001) integral field spectroscopy. The stellar and gas 
kinematics has been published in Ganda et al. (2006), while the ab- 
sorption line maps have been analyzed in Ganda et al. (2007). Our 
sample of late-type galaxies can be easily compared to the early- 
type spirals of the SAURON Survey (de Zeeuw et al. 2002), pre- 
sented in Falcon-Barroso et al. (2006) and Peletier et al. (2007). The 
fact that kinematic parameters are available for this sample makes it 
possible to study the position of late-type spirals on the M^suige - o" 
relations, establishing in this way whether these bulges are similar 
in structure and formation to bulges of early-type spirals. 

Given the renewed interest in the extinction in spiral galaxies 
(Driver et al. 2007, Graham & Worley 2008, hereinafter GW) we 
decided to use the unique property of this sample that the amount 
of extinction can be calculated by combining absorption line maps, 
which are basically independent of extinction, with broad band 
colour maps, such as the Mg b maps, which are affected by it. To do 
this, we determined high resolution V ~ H maps for the 10 galax- 
ies for which HST-F606W {V) images were available, in addition 
to F160W (H). The combination of colour and line strength gives 
the colour excess YLv-h in a model-independent way, which can 
easily be converted to Av - This method is in principle very power- 
ful, but has not been applied much in the literature, because of the 
lack of well-calibrated line strength data. 

The paper is structured as follows. Section |2] describes the 
sample selection; Section|3]the photometric profiles and the meth- 
ods applied for their extraction; Section|4]the actual bulge-disk de- 
composition; Section [5] compares the photometry in the NIR with 
the optical V~ band, and presents extinction maps; Section [6] in- 
vestigates the correlations between the structural parameters and 
other galaxy properties; Section [7] addresses the frequent presence 
of inner additional components, and, finally, Section[8]summarizes 
the main results. 



2 THE SAMPLE 

The sample of galaxies on which we perform our bulge-disk 
decomposition is the same for which we presented the two- 
dimensional kinematical maps from SAURO N integral-field spec- 
troscopy in the paper by Ganda et al. and the line-strengths 
maps in the paper by Ganda et al. (20()7|). 

The galaxies wer e optically selected (Bt < 12.5, according 
to the values given in lde Vaucouleurs et al.|[l99ll hereafter RC3) 
with HST imaging available from WFPC2 and/or NICMOS. Their 
morphological type ranges between Sb and Sd, following the clas- 
sification given in NED (from the RC3). Galaxies in close interac- 
tion and Seyferts were discarded. The resulting sample contains 18 
nearby galaxies, whose prope rties are listed and illustrated respec- 
tively in Table 1 and Fig. 1 in lGanda et al. I dioo^ . 



3 PHOTOMETRIC PROFILES 
3.1 Datasets used 

We constructed surface brightness profiles for all galaxies combin- 
ing datasets with different spatial resolution and field of view. We 
started from the H~ band images from the Two-Micron All Sky 
Survey (hereafter, 2MASS). Since these data do not have a suffi- 
cient spatial resolution for our purposes (2-3"; the pixel scale of 
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the 2MASS images retrieved from the archive is 1"), we comple- 
mented the photometric profiles extracted from the 2MASS images 
with profiles extracted from NIR HST images (F160W), available 
for 11 out of 18 cases (marked with a single or double dagger in 
Table [T](. To improve the spatial resolution at small radii for the 
remaining galaxies we used optical HST images (F814W, corre- 
sponding to the /— band). Since, on the other hand, 2MASS is 
rather shallow (reliable profiles can be extracted out to ~ 5 - 8 
kpc for most of our galaxies), we needed an additional dataset for 
the outer parts, so that the disk geometry could be determined un- 
ambiguously, and an accurate sky level could be determined. For 
this reason, we decided to complement the profiles also with pro- 
files extracted from infrarecQ images from the Digitized Sky Sur- 
vey (hereafter, DSS), which are very deep, and therefore radially 
very extended. This also allows us a more accurate determination 
of the sky background. The DSS, 2MASS and NICMOS images 
downloaded from the archives were fully reduced, so we did not 
apply any further processing; as for the WFPC2-F814W images, 
we first created a mosaic from the four chips in order to maxi- 
mize the field of view, using the IRAF task WMOSAIC in the STS- 
DAS.HST.CALIB.WFPC package, rotated the resulting image to ori- 
ent it to North up-East left and performed cosmic rays removal, 
using the IRAF tasks CRREJ and LACOS.IM, created by Pieter van 
Dokkum and available via his websitfl 



3.2 Ellipse fitting 

As a first step, we extracted the photometric profiles from 
the DSS images using the IRAF task ELLIPSE in the STS- 
DAS.ANALYSIS.ISOPHOTE package, which fits elliptical isophotes 
to galaxy images, impl ementing the method initially described by 
Kent (1983, 1984) and ' jedrzeiewskil ( Il987l) . We masked out bad 
pixels and foreground stars, but did not exclude dust lanes and star 
forming regions from the fitting. We first fitted ellipses to the im- 
ages with the centre, position angle and ellipticity left as free pa- 
rameters. In this way we obtained profiles of the centre coordinates; 
from these profiles we established the position of the centres, kept 
fixed in the following steps. In case of ambiguity, the centre was de- 
termined using isophotes at intermediate radii, since the very cen- 
tral ones can be affected by dust obscuration, and, in the particular 
case of the DSS images, by the seeing (or even by saturation of the 
images, in the worst cases). 

We then fitted again ellipses to the images with the centre fixed 
at the chosen position and ellipticity and position angle free. From 
this second step we determined single values for the geometric pa- 
rameters. In many cases in the outer parts of the galaxies, outside 
the nuclear region, the position angle and ellipticity converge to 
more or the less constant values. In other cases the situation is 
more complicated, especially for galaxies that are almost round, 
and there we picked the values that looked more reasonable at a 
visual inspection of the shape of the isophotes in the images. 

Table [T] lists the position angles and ellipticities of our sam- 
ple galaxies, together with some other basic data, and shows that in 
most cases the chosen values are very close to the values tabulated 
in the literature (mainly in the RC3); the cases which present the 
biggest discrepancies are actually galaxies almost face-on and/or 



1 We used POSS-H Near-IR IVN+RG9 images, obtained from 
lhffpi77archive.stsci.edu/cgi-bin/dss_form/ the selected filter is close to 
the Z~ band. 

^ http://www.astro.yale.edu/dokkum/lacosmic/ 
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Table 1. Column (1): NGC identifier; column (2): morphological type 
(RC3); column (3): absolute blue magnitude, from HyperLeda, computed 
using distance in Mpc given in column (4); column (5): scale in pc/arcsec; 
columns (6) and (7): adopted position angle (PA) and ellipticity (e) for the 
measurement of the photometric profiles (N-^E); columns (8) and (9): po- 
sition angle and ellipticity values tabulated in the literature (PAm and e;it); 
we refer to the RC3 values, with the exceptions of NGC 1042, 3346, 4254, 
4775, 5668, where the angle is taken from Grosb0l 1 1985). The galaxies 
marked with a single or double dagger are those for which we used NIC- 
MOS imaging from the NIC2 or NIC3 cameras, respectively. For the re- 
maining galaxies WFPC2-F814W images are available. 



very round, for which the errors on the geometric parameters are 
in any case big. Once we had determined the geometric parame- 
ters, using ELLIPSE again we measured the intensities at each ra- 
dius with centre coordinates, position angle and ellipticities fixed 
for all ellipses. From the resulting photometric profiles we sub- 
tracted then a value for the sky background estimated on the outer 
parts of the images. Following this, we fitted ellipses to the 2MASS 
images, with all the geometric parameters (centre coordinates, po- 
sition angle and ellipticity) free and determined the centre coordi- 
nates in the same way as done for the DSS images. And after this 
we measured the intensity profiles, fixing the centre coordinates 
to the values just determined and the position angle and elliptic- 
ity to the values determined on the basis of the DSS profiles and 
subtracted an estimated value for the sky brightness. An analogous 
procedure was repeated for the NICMOS images, when available, 
and for the WFPC2-F814W images, in the cases where NICMOS 
was missing (see Table[T). 



3.3 Combining the profiles 

At the end of this procedure, for each galaxy we combined the three 
photometric profiles, to get a single profile with the maximum ra- 
dial extension and spatial resolution allowed by the data. To do this, 
we first combined the profiles extracted from the DSS and 2MASS 
images, aligning them at the 2MASS level, and transformed the re- 
sulting 'ground-based' profile to an absolute magnitude scale using 
the zero point calibration (the keyword MAGZP) of the 2MASS 
images. For the 1 1 galaxies with available NICMOS imaging, we 
first calibrated independently the NICMOS profile to the absolute 
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H (Vega) magnitude scale using the F160W zero points of 21.826 
(plus a term taking into account the pixel scale)for the 10 galaxies 
with NIC2 data and 21.566 for NGC 5585, observed with NIC3, 
and then joined them with the calibrated 'ground-based' profiles; 
the quality of the matching of the two profiles gave us an indica- 
tion of the quality of our estimation of the sky background in the 
HST image, which was then changed iteratively until a satisfac- 
tory result was obtained. For the remaining galaxies, we imposed 
as a photometric zero point to the HST profiles (F814W) the mean 
offset between the aperture magnitudes measured on the 2MASS 
and HST images in the radial range 7-12", which depends on an 
estimate of the sky background level for the HST image; we then 
joined the HST and the calibrated 'ground-based' profiles, modi- 
fying our sky estimate for the HST image until a good match was 
obtained; given the fact that the field of view of WFPC2 is sig- 
nificantly larger than the one of NICMOS, this operation required 
much less careful fine-tuning than for the NICMOS profiles. At the 
end, calibrated global H— band profiles were obtained. We tested 
whether the use of WFPC2-F814W images in the cases where NIC- 
MOS is not available introduces significant additional errors: for 
several of the galaxies with available NICMOS imaging, we re- 
trieved the F814W images and treated them in the way described 
above; the result is that the comparison between the final global 
profiles obtained this way and the ones obtained using NICMOS 
is acceptable, with differences below 0.1 mag outside the inner ~ 
l'.'5 and below 0.05 mag outside the inner ~ 4'.'5. This means that 
we can be rather confident in using the F8 14W images for the inner 
regions, and that I — H colour gradients are not very important. 



4 BULGE-DISK DECOMPOSITION 

Fitting a parametric bulge and disk can be done directly on the 2- 
dimensional images, or on one-dimensional profiles, that have been 
obtained by averaging azimuthally in ellipses or circles. Although 
one-dimensional profiles suffer from loss of information due to pro- 
jection, and produce larger degeneracies between bulge and disk 
parameters l lde Jond ( fT996ah), two-dimensional fits are less robust, 
and degeneracies remain TMollenhoff & Heidt"200l'). On the ba- 
sis of their simulations MacArthur et al. ( 2003) conclude that a 
full two-dimensional technique does not provide a significant im- 
provement compared to one-dimensional methods in recovering the 
axysimmetric structural parameters. Based on this, we have chosen 
an approach that allows us not to have to assume the same shape 
of bulge and disk, but for which we do not have to apply a full 
two-dimensional fitting method, as we will explain in the follow- 
ing sections. Alternatively, we could have used a non-parametric 
fitting method, such as the method of Kent (1984), where two com- 
ponents are found with different axis ratios, but have not done this, 
because of the degeneracies that arise when choosing the axis ratios 
of these components. 

4.1 Decomposing late-type spirals 

We performed a non-simultaneous bulge-disk decomposition for 
our 18 galaxies, c lose ly following the method de s cribed by 
iNoordermeeil ( |2006|) and Noordermeer & van der Hulst ('2006") and 
previously introduced by Palunas & Williams (2000). In a few 
words, we fit an exponential law to the disk, marking the fitting 
range via a visual inspection of the profiles, choosing the part of 
the profile that is exponential without doubt; we build a model im- 
age of the disk and subtract it from the original images, obtaining 



'bulge images' from which we extract the 'bulge profiles' that we 
then fit with a Sersic r^''" law. 

4.1.1 disk fitting 

The first step in our bulge-disk decomposition is the estimation 
of the disk parameters by means of profile fitting. On the com- 
bined global profiles, measured as explained in Section [3] using 
a least-square algorithm we fitted an exponential light distribution 
(in mag arcsec"^): 

=Mo,<i + 1.0857 (^0 (1) 

where /.io,d and h are the central surface brightness and scale 
length of the disk. The actual fit was performed on the profiles 
in absolute magnitudes and on radial ranges chosen by means of 
a visual inspection, in order to avoid the contamination from the 
central component (the 'bulge'). The choice of a fitting range be- 
came a sort of 'compromise' in a few cases where the shape of the 
profile would have su ggested the presence of a double disk (see 
IPohlen & Truiilldl2006 ). Fig. [T] shows the global galaxy profiles, 
together with the best-fitting exponentials, and the residuals from 
the fit. The radial ranges of the fits are indicated in the residual plot. 
These figures show that in general the fits are satisfying. NGC 4102 
represents the case where the fit is worst; this might be due to con- 
tamination from the big bar, or from the coexistence of an inner 
and an outer disk (see also Pohlen & Trujillo 2006), to which we 
impose a single-exponential fit. 

The fit parameters are listed in Table [2] No correction for in- 
clination or galactic extinction is applied. The uncertainties of these 
parameters can be calculated in various ways. Often uncertainties 
in the least-squares fit are quoted. These are often very small, since 
galactic disks cane be fitted very well with an exponential surface 
brightness profile, if a limited radial fitting range is used (de Vau- 
couleurs 1959, Freeman 1970, MacArthur et al. 2003, Pohlen & 
Trujillo 2006). By changing the radial range, however, disk param- 
eters can vary considerably. Therefore, a more realistic estimate of 
the uncertainties in surface brightness and scale length is obtained 
if one studies the variation of these parameters as a result of chang- 
ing the fitting range. To do this, we varied both the inner and outer 
radial limit by steps of 5%, from 70% to 130% of the range quoted 
in Table 2. The final parameters disk are the average parameters 
obtained in this way, and the uncertainty the RMS scatter. 

4.1.2 Bulge fitting 

In all cases, in the innermost regions the galaxy's light clearly ex- 
ceeds the exponential fit to the disk. According to the photomet- 
ric definition (see the Introduction), the region of the galaxy corre- 
sponding to this light excess is the bulge. We will adopt this defini- 
tion, but without implications on the kinematics and/or populations: 
without implying that it has to be a hot and old component, similar 
to a small elliptical. 

To quantify the region where the light from the bulge is sig- 
nificant, we established a bulge radial extension, defining as bulge 
the region within which the intensity ratio between the fitted expo- 
nential disk and the total intensity is lower than 0.7 (i.e. the bulge 
contributes to 30 percent or more of the total light, see Fig. [T] and 
Table|2}. Note that these figures can be affected by the presence of 
bars or other asymmetric features not belonging to the exponential 
disk. 

On the basis of the fitted disk parameters, using the standard 



The nature of late-type spiral galaxies: structural parameters, colour profiles, and dust extinction 5 



12.0 
16.0 
20.0 

24.0 
0.4 
0.0 

-0.4 
2.0 
1.0 



NGC488 
SA(r)b 



NGC2964 
SAB(r)bc 



NGC4254 
SA(s)c 




NGC772 
SA(s)b 




NGC3423 
SA(s)cd 




NGC4775 
SA{s)d 



: NGC864 : 
1 SAB(rs)c : 


NGC3949 - 
SA(s)bc : 


NGC5585 : 
SAB(s)d : 


^^*t^.^,.t| : 


- #*^hV+-^+^+++^+- 

J \ w- - 




^^**H^-iH^±+_+^+tt : 


^ r- : • ! + + + ^ + 


|,^^Hi+t++++++ti± + + + ^ + - 



"Eb 

n ■ 

s 



fc: 
o 



NGC5668 
SA(s)d 




50 100 150 200 250 300 
Semi-major Axis [arcsec] 



50 100 150 

Semi-major Axis [arcsec] 



50 100 150 200 

Semi-major Axis [arcsec] 



Figure 1. Global combined photometric profiles and exponential fits to the disks. Filled dots show the global combined photometric profiles, along with the 
best-fitting exponential (solid line). The units are arcsec along the semi-major axis of the ellipses on the horizontal axis and absolute H magnitudes arcsec^^ 
along the vertical. Below each profile the residuals from the fit are shown, in the same units. Below this, the disk to total ratio of the intensity of the fitted 
exponential disk and the total profile intensity is shown; the level D jT = 0.7 is indicated (dashed line), defining the 'bulge region', as explained in the text. 
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IRAF tasks (BMODEL in the STSDAS. ANALYSIS. ISOPHOTE pack- 
age) we built a model image of the disk and subtracted it from the 
original 2MASS images. We will refer to the residual images as the 
2MASS 'bulge image'. From these images we extracted the bulge 
photometric profiles and parameters, by applying a procedure sim- 
ilar to the one used to measure the 'total' profiles as explained in 
Section[3] we fitted ellipses to the bulge images in two steps, keep- 
ing in all cases the centre fixed to the value already determined 
for the original 2MASS images and the position angle fixed to the 
value determined from the DSS images in Section[3] At first, we left 
the ellipticity free; in a second step we fixed also the ellipticity and 
extracted the intensity profiles from the 2MASS bulge images, suc- 
cessively converting them to an absolute magnitude scale, applying 
the offset given by the keyword MAGZP in the 2MASS headers. 
As ellipticity for the bulge we adopted a representative value esti- 
mated within the bulge region, on the basis of the ellipticity profiles 
extracted from the 2MASS and HST images. We notice that bulge 
and disk ellipticities can be different. Fig. [2] presents a compari- 
son between the adopted values for the ellipticities of disks and 
bulges, showing that for most of the galaxies the inner component 
is rounder than the disk; but in some cases (NGC 3346, 4487, 4775, 
5668, among the latest-type objects in our sample), the bulge is in- 
stead flatter than the disk, as shown also bv lFathi & Peletie 3 ^20031) 
on the basis of the bulge-disk decomposition of 70 disk galaxies 
spanning a range in type between SO and Sm. Only NGC 3346 lies 
more than \a above the 1:1 relation in Fig. (2] 

Once the bulge and disk profiles have been retrieved, we need 
to determine the structural parameters by fitting a model profile. We 
follow the approach most often used in the current literature (see 
Graham & Worley 2008 for a comprehensive review) by fitting an 
exponential disk and a Sersic r^^" bulge. 

When fitting the profiles, seeing effects are particularly rele- 
vant when the effective radius is small and/or when the ratio be- 
tween the effective radius and the FWHM of the seeing is small 
dOra ham 2001). In fact, for small bulges poor seeing could smear 
the images so that bulges intrinsically described by n > 1 Sersic 
profiles could appear like n < 1 profiles. Therefore, for a reliable 
treatment of the seeing effects, we performed the fit on the bulge 
region by convolving the Sersic profiles with a Gaussian point- 
spread function, as explained by Graham (2001), and then fitting 
the seeing-convolved Sersic profiles using a standard non-linear 
least-squares algorithm. The values of the FWHM of seeing in the 
2MASS images were retrieved from the header keyword SEESH 
via the relation FWHM(") = 3. 13 x SEESH - 0.46 as explained on 
the 2MASS webpagtH. In some cases (e.g. NGC 864, 2964, 3346, 
3949, 4775 and 5678) the bulge region is so small that it is com- 
parable in size with the FWHM of seeing in the 2MASS images 
(less than 3 x FWHM of the seeing), so that the applied proce- 
dure is not meaningful and possibly leads to imprecise fitting val- 
ues. In order to face this problem, to gain spatial resolution in the 
inner parts, and to take advantage of data non contaminated by see- 
ing effects, we perform the bulge fitting on bulge profiles extracted 
from the HST images (NICMOS or F814W, depending on avail- 
ability). These profiles were obtained by subtracting from the HST 
images an exponential disk, modeled on the basis of the fit param- 
eters obtained as described in Section 14.1.11 and then fitting to the 
residual images ellipses with position angle and ellipticity fixed to 
the 'disk position angle' and the 'bulge ellipticity', and centre co- 
ordinates fixed as well. Some of these 'HST bulge profiles' suffer 
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Figure 2. Bulge against disk ellipticities; each filled dot represents a galaxy, 
with the NGC number indicated close to the symbol. The disk ellipticity is 
determined from the outer parts of the galaxies, the bulge ellipticity from 
the inner regions, as explained in the text. The dotted line represents the 
1:1 relation. The overplotted error bars are obtained from the scatter of the 
points in the profiles. 



from the complementary problem of limited spatial extension: in 
several cases they do not cover the full bulge region. In those cases 
(NGC 488, 628, 772, 4030, 4254, 4775, 5668) we joined them with 
the bulge profiles previously obtained from the 2MASS images. 
On these profiles we fitted a pure Sersic law, excluding the very 
innermost parts (0.5-1"), which in several cases host nuclear star 
clusters. In Fig.[3]we show the HST or joined HST -I- 2MASS bulge 
profiles, along with the best-fitting Sersic profile, and the residuals 
from the fit. 

The fit parameters for all galaxies, for both the disk and the 
bulge components, are listed in Table|2] Here we report the disk pa- 
rameters from the exponential fit as in Section |4. 1.11 for the bulge 
we quote the parameters from the seeing-convolved Sersic fit to the 
2MASS bulge profiles and those from the Sersic fit to the HST or 
combined 2MASS-HST bulge profiles. No correction for inclina- 
tion or galactic extinction is applied to the surface brightness (and 
magnitudes) values reported in the table. The uncertainties of the 
bulge parameters have been calculated in the same way as for the 
disk, i.e. by varying the fitting range by 30% for both the inner and 
outer boundary. Fig.|4]shows the comparison between the two sets 
of parameters for the bulge; the effective surface brightness and ra- 
dius are generally rather similar, while the Sersic parameter n tends 
to be larger for the fits on the HST (or combined HST -I- 2MASS) 
profiles. In the following, for the bulge parameters we will refer to 
the output of the fit on these latter profiles. From the parameters we 
computed the total bulge and disk luminosities (see Graham 2001), 
and converted them to a bsolute magnitudes M^ and M;, using the 
distance moduli used by iGanda et alj ( l2006h . These numbers, to- 
gether with the bulge to disk ratio, are also given in Table|2] 
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Figure 3. Sersic fits to tlie bulge profiles. In the top part of each panel we plot the HST or combined HST + 2MASS bulge profiles (filled dots), in mag 
arcsec" ; the NICMOS images are used in all available cases, the F814W-WFPC2 in the others (see TableQ; overplotted are the best-fitting Sersic profiles 
(solid line); The low part of each panel shows the residuals from the fit (in mag arcsec"^), as well as the fitting range (vertical lines). Note that the abscis is 
given in logarithmic units, contrary to Fig.[T] 
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Figure 4. Comparison between the fitted pai'ameters for a seeing-convolved Sersic fit performed on the 2MASS bulge profiles and a pure-Sersic fit on the HST 
or joined HST + 2MASS bulge profiles, excluding the innermost region: from left to right, effective surface brightness, effective radius, n. The dotted line in 
each panel represents the 1 : 1 relation. 
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Table 2. Column (1): NGC identifier; columns (2) and (3): disk fitting range, in arcsec; columns (4) and (6): respectively, central surface brightness of the disk 
in mag arcsec""^ and scale length of the disk in arcsec, from the fits performed on the photometric profiles in Section l4.1.1l column (5) and (7) indicate their 
uncertainties; columns (8) and (9): bulge fitting range, in arcsec; columns (10), (11) and (12): effective surface brightness in mag arcsec^^, effective radius in 
arcsec and Sersic parameter (n) of the bulge from the fit of a seeing-convolved Sersic profile to the bulge profiles extracted from the 2MASS images; columns 
(13), (15) and (17): bulge parameters (effective surface brightness, effective radius and Sersic parameter (n) from the fit of a Sersic profile to the HST or joined 
HST -I- 2MASS bulge profiles, with uncertainties in columns (14), (16) and (18); columns (19) and (20): H— band absolute magnitudes of the disk and bulge 
components, and the uncertainty in the latter (21). The uncertainties quoted here are the variations due to changing the fitting ranges. 



5 OPTICAL PHOTOMETRY AND OPTICAL - NIR 
COLOURS 

As we mentioned in the Introduction and as we witnessed through- 
out this work, late-type galaxies are known to be dusty objects, with 
dust lanes an d structures extending often all the way to the centre 
jFrogell 19851. Zaritsky, Rix & Rieke 1993, Martini et al. 2003). For 
example we just mentioned in the previous section that dust makes 
it sometimes difficult to investigate the properties of tiny inner com- 
ponents. This is one of the reasons for which we decided to work 
in the near-infrared. In order to investigate the amount and distri- 
bution of extinction in the centers of these objects, we decided to 



analyse also optical images. Here, the fact that we have absorption 
line strength maps determined with SAURON offers a large advan- 
tage, when compared to previous work. Comparing optical-near in- 
frared colours, affected by extinction, with absorption line index 
maps, which are almost not affected, we can obtain a reasonably 
accurate guess of the amount of extinction in the central regions of 
a representative sample of late type spirals. There are 13 galaxies 
for which F606W HST images are available in the HST archive. 
For 10 of them, also NICMOS F160W images are available. We 
decided to analyse these 10 galaxies. 
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5.1 Optical - Near-infrared colour profiles 

We analyzed the optical images in the same way as the NIR ones, 
with the exception that here we used the POSS-II Red IIIaF + 
RG610 images in the outer parts: using the IRAF task ELLIPSE 
we first extracted from the DSS images photometric profiles with 
centres and geometric parameters free and determined the centres 
coordinates; we then fixed the centres and extracted again profiles 
with free ellipticities and position angles. In many cases in the outer 
parts the geometric parameters converge to a constant value; in the 
majority of cases this value is comparable with what found in the 
NIR: we did not detect any significant difference in the shape of 
isophotes between optical and near-infrared. Therefore, we fixed 
the position angle and ellipticity to the values determined on the 
NIR images (see Table [TJ and measured the intensities along the 
ellipses. We then subtracted the sky background, estimated by av- 
eraging the mean values of intensity in several boxes placed in the 
outer parts of the images, and converted the profiles to a (uncali- 
brated) magnitude scale. For the HST F606W images, we started by 
creating mosaic images of the four chips, using the IRAF task WMO- 
SAIC, rotated them to bring them to the standard orientation with 
North up-East left and extracted from them intensity profiles with 
centre, position angle and ellipticity fixed (fixed to the NIR values); 
we then subtracted an estimate for the sky background and con- 
verted the profiles to a (still uncalibrated) magnitude scale. Then 
we combined HST and DSS profiles into single profiles, aligning 
them at the level of the HST profiles; when a good matching of the 
two was not achieved, we modified the estimated sky brightness in 
the HST image, until the matching was satisfactory. The last step 
is the photometric calibration: converting the profiles to an abso- 
lute magnitude scale. We calibrated our profiles to the V— band of 
th e UBVRI Johnson -Cousin system using the equations provided 
bv lHoItzmanl ( ll995h . that were derived for single-chip images. To 
be able to apply those transformations, for each galaxy we focused 
on the single-chip image containing the galaxy centre (PCI for all 
galaxies excluding NGC 1042 and 3949, centred instead on WF3) 
and extracted the single-chip profiles. The sky subtraction is in this 
case more difficult, since the field of view is significantly smaller: 
we proceeded in an iterative way, changing the estimated value 
of the sky background until the sky-subtracted single-chip profile 
matched the shape of the sky-subtracted profile retrieved from the 
mosaic image . We then calibrat ed the single-chip profiles applying 
Equation 9 in iHoltzmad fT99^ . using the coefficients listed in his 
Table 10: 

fiv = -2.5 X log{counts) + 2.5 x log{EX PTI M E) + 

+ 5 X \og{SCALE) + 0.254 x (V - /) + 0.012 x (2) 

x{V - if + 22.093 + 2.5 x \og{GR) + 0.1 

where 'counts' is the number of counts in linear units, 
'EXPTIME' the integration time on the image, 'SCALE' 
the pixel scale (arcsec pixen^, ~ 0.046 for the galaxies centred 
on PCI and ^ 0.09944 for the two centred on WF3), V - I 
the colour in the Johnson-Cousins UBVRI system , 22.093 the 
photometric zero point tabulated bv lHoltzmanI il995l) , GR equals 
2.003 for frames where a gain of 14 was used , and 1 otherwise , 
and 0.1 is added, following the prescription of iHoltzmanI ( Il995h . 
in order to correct for infinite aperture (see also the WFPC2 Data 
Handbook). For the (V — I) colours of our galaxies we used 
values from HyperLeda and, when not available, an average of 
the values for galaxies of the same Hubble type. Since the colour 
term is relatively small, there is no need for extremely accurate 
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Table 3. Values for the colour V — I adopted in the calibration of the optical 
images from HST. 

values here. The adopted values are listed in Table |4] The last step 
was to find and apply the vertical offset necessary to align the 
HST(mosaic) + DSS joined profiles to the calibrated single-chip 
profiles. 

Four of the galaxies for which we perform our optical 
analysis (NGC 104 2, 3423, 4030 ari d 4102 ) belong also to the 
sample studied by IPohlen & Truiilld ( l2006h . who present SDSS 
g' and r' profiles (roughly corresponding to B and R bands) for 
90 almost face-on late-type galaxies. The authors kindly provided 
us their profiles for the galaxies in common. In their profiles we 
recognize all the features that we see in ours, at the same radii, 
out to ~ 150 - 200" (depending on the galaxy): the profiles are 
completely consistent, with the caveat that they are measured in 
different bands. 

From the thus obtained V and H profiles we determined ra- 
dial V — H colour profiles, that we present in Fig.|6] Here we only 
present the inner 20", since real near-infrared imaging is only avail- 
able in this region. Further out, the 2MASS images are not deep 
enough. We find that the profiles have large dispersions in colour, 
and a wide variety in slope. Some galaxies become redder, some 
bluer as a function of radius. When comparing them to elliptical 
galaxies (e.g. Peletier, Valentijn & Jameson 1990) the profiles are 
much less smooth, and the gradients much larger. The colour maps 
are very instructive in telling us what is happening. Fig.|7ji and|7j5 
show the calibrated V — H colour maps of the inner 20" x 20" 
of the galaxies. They have been displayed on the same scale, to 
make it easier to make the comparison from galaxy to galaxy, and 
to compare the internal gradients with the colour differences from 
galaxy to galaxy. The colour maps show dust lanes in all galaxies 
except for NGC 5585, spiral arms with younger stars, inner spirals, 
nuclear clusters, etc. In most cases it is not easy to find out where 
the galaxy center is, when just looking at the colour map. 



5.2 Extinction in the centers of late-type spiral galaxies 

From aV — H map alone it is very difficult to determine the amount 
of extinction Ev-h- The problems is that the intrinsic colour of the 
stars is unknown, so that the uncertainty in the colour of the stars 
is reflected in the uncertainty in the amount of extinction. Having 
two colours does not help very much, since the effects of redden- 
ing in colour-colour diagrams is almost parallel to the effect due to 
changing metallicity or age (e.g. Kuchinski et al. 1998), and since 
the distribution of the dust changes the relative extinction in the 
different bands. Determining the amount of extinction in a specific 
galaxy is therefore done assuming a certain colour for the stellar 
populations (see e.g. Knapen et al. 1995), using the Balmer decre- 
ment for ionised gas (which usually gives the extinction on the li- 
neof sight to an HII region, but not in the whole galaxy), or using 
statistical methods. For example, from the distribution of colour 
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Figure 5. V - H colour profiles for the 10 galaxies with available optical and NIR HST images 
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Figure (t.V — H colour vs. Mg b line strength in central apertures of 2.4" for the 10 galaxies with NICMOS F160W and F606W HST data (black filled 
symbols). In red (asteriscs) elliptical and SO galaxies are shown from Kuntschner et al. (2006), with Y — K measurements from Frogel et al. (1978). Lines 
represent SSP models of Bruzual & Chariot (2003) for various metallicities: green (long dashed-dotted) indicates Z=0.05, blue, full lines indicate Z=0.02, red 
dotted line Z=0.008 and magenta short-dashed-dotted lines Z=0.004. 



as a function of inclination one can derive the amount of extinc- 
tion (e.g. TuUy et al. 1998, Giovanelli et al. 1995, Peletier et al. 
1995). The SAURON dataset offers the nice advantage that absorp- 
tion line indices are available in the inner 40" x 30". Since to first 
order the Mg b index in a galaxy gives the same information as the 
Y — H colour (which is the metallicity for an old galaxy, see Fig. 
[6](, and since Mg b is affected very little by extinction (MacArthur 
2005), one can use the V H and Mg b together to determine 
the amount of extinction at every position where both are avail- 



able. Here we use the central Mg b measurements of Ganda et al. 
(2007) to calculate the extinction in the central aperture with di- 
ameter 2.4". In Fig. [6] these central values are shown. The figure 
shows that the models with various metallicities almost fall on top 
of each other. This means that one can measure the extinction by 
measuring the distance in V ~ _ff to a line of models with a given 
reference metallicity. Since most of the galaxies have luminosities 
similar to the Milky Way, we take the solar models as a reference. 
For NGC 5585, which is fainter, it is probably more appropriate to 
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Figure 7. (a): Extinction distribution for the 10 galaxies with space-based photometry in F606W and F160W. Left: surface brightness distribution in H. Middle: 
calibrated V-H colour maps. Right: extinction maps from V — H and Mg b. The ranges of the plots are shown in the lower right of each column of plots. 




Figure 7. (b): Extinction distribution for the 10 galaxies with space-based photometry in F606W and F160W. Left: surface brightness distribution in H. Middle: 
calibrated V-H colour maps. Right: extinction maps from V — H and Mg b. 
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3.16 


0.00 


0.00 


4030 


3.20 


0.47 


0.57 


628 


2.92 


0.33 


0.40 


4102 


3.76 


1.44 


1.75 


864 


2.57 


0.35 


0.42 


4254 


3.32 


0.74 


0.90 


2964 


3.27 


0.90 


1.09 


5585 


2.24 


0.22 


0.27 


3949 


2.52 


0.14 


0.17 


5678 


3.60 


1.05 


1.27 



Table 4. V — and central extinction in an aperture of diameter 2.4" 
derived from V — H and Mg b (see Fig. [6). 



take a lower metallicity. The uncertainty in the amount of extinc- 
tion is then given by the uncertainty in the V ~ H of the stellar 
populations in the galaxy that make up the Mg b. From the dis- 
persion between the SSP models we estimate that this uncertainty 
is about 0.2 mag or less, even if the stellar populations consist of a 
mix of ages, since more complicated models are always linear com- 
binations of SSPs. The obtained Ev~h values are converted to Av 
using the Galactic extinction law (Rieke & Lebofsky 1985), divid- 
ing them by 0.825. Here the assumption is made that the galaxies 
are not optically thick in H. Given the extinction values that are ob- 
tained, this assumption is easily satisfied. Table 4 shows the extinc- 
tion values we derived. We find an average Ev^h of 0.56 mag, or 
a lower limit to the average extinction Av=0.68 mag. In Fig.Qthe 
central extinction maps, determined in this way, at the SAURON 
resolution, are shown. 

Do we learn anything about the stellar populations from the 

V — H colour? Figure [6] shows that we don't. For some galaxies 

V — H is redder than for the most massive giant ellipticals. Since 
the line indices of all objects are weaker than those of massive gi- 
ant ellipticals, we have a strong indication that it is dust extinc- 
tion, and not large numbers of very red stars, that causes the red 

V — H colour. This argument is made stronger by the fact that for 
highly inclined galaxies optical-infrared colours in regions outside 
the dustlane always are bluer than those of giant ellipticals (Peletier 
& de Grijs 1998). The issue of the very red optical-infrared colours 
was already discussed in 1985 by Frogel, who found galaxies with 

V — K colours that were considerably redder than those of the 
brightest giant ellipticals. He found also that half of his sample of 
Sc galaxies had U — V colours that were bluer than the giant el- 
lipticals. Although he could not uniquely determine the extinction, 
he used the following experiment. If all galaxies had an intrinsic 
stellar population colour of J — H=0.65, he obtained extinction 
values between and 1.64 mag in Av. The optical-infrared colours 
would not be redder any more than those of giant ellipticals, and 
also reasonable U ~ V colours would be obtained. Independently, 
Turnrose (1976) had derived internal extinction by comparing sim- 
ple stellar population models with observed optical stellar energy 
distributions, and had found similar extinction values. We have one 
galaxy in common with Frogel (1985), NGC 628, for which V — H 
agrees within 0. 1 mag. 

Since these are all close to face-on galaxies, the measured ex- 
tinctions are not very different from what one would see if they 
were face-on. The extinction correction to face-on is given by, 

Afiv = 7 log{a/b) 

where a/b is the axis ratio of the galaxy, and 7 a transparency index 
(Boselli & Gavazzi 1994). For this sample log{a/b) ranges from 
to 0.26, while values for 7 should be between 0.4 and 0.7 (Boselli 
& Gavazzi 1994), so that the face-on corrected Av is between and 
0.2 mag less than the value given in Table|4l sometimes still around 
or above 1 mag. This result is compatible with TuUy et al. (1998), 



Peletier et al. (1995) and Giovanelli et al. (1995), who used a sta- 
tistical method, deriving the extinction from the relation of colour, 
or surface brightness, as a function of inclination. 

We now look at extinction as a function of morphological type. 
This is shown in Fig. [8] The relation is not very clear, but it can be 
understood when one looks at the outliers. The galaxy of morpho- 
logical type 7 (Sd, NGC 5585) is by far the bluest, and the small- 
est and has little extinction. This is in agreement with Tully et al. 
(1998) who find that galaxies with Mb =-18 have almost no extinc- 
tion. At the other end, at morphological type 3=Sb, we find a galaxy 
with very little extinction, NGC 488, and two others with a lot of 
extinction. This can be understood by looking at Fig. 12 of Tern- 
drup et al. (1994), where one can see that galaxies with type 3-5 
have much more extinction than galaxies of earlier types. Peletier 
& Balcells (1996) confirm this, showing also that the extinction in 
galaxies of type 7 and larger is again much smaller than in galax- 
ies of type 3-6. It looks as if NGC 488 is more similar to an Sab 
galaxy than an Sb in its inner regions. One should note that the Av 
values that we have derived here are still lower limits. For example, 
by adding very young stellar populations mixed with dust, one can 
still make the color of a galaxy bluer. For that reason, the conclu- 
sion that the center of NGC 3949 is almost dust free is most likely 
wrong looking at the structure in the colour map. Our derived ex- 
tinction values for these late-type spirals are consistent with Tuffs 
et al. (2004), who use the models of Popescu et al. (2000) to de- 
termine the attenuation of dust in spirals using a thin and a normal 
disk with extinction, and a bulge without dust. We should note that 
this models has a very large optical depth (TB,cen ~ 3.8) so that, 
according to the model, we basically only see the bulge light in 
front of the disk, as a result of which the integrated reddening is 
not so high, and the colours reasonable. With these data we are not 
able to verify whether this extinction model is realistic for late-type 
spirals. 

Does the finding that many of our late-type galaxies have a 
considerable amount of extinction in their centre have other con- 
sequences? Martini et al. (2003), and Carollo et al. (2002) also 
show optical-infrared maps of similar galaxies, some of which are 
in common with us. They however do not show any calibrated color 
maps. To test the effects of extinction on the structural parameters, 
we also calculated the same parameters as presented in Section 5 
in the V-band. In general, the differences were small, something 
which could have been expected from the fact that for these galax- 
ies the colour profiles are generally reasonably flat. Since there are 
many galaxies with large colour gradients (see e.g. Martini et al. 
2003 and Carollo et al. 2002) we think it is preferable to determine 
these parameters in the NIR only. 

Peletier & Balcells (1996) show that colours of bulges and 
disks of spirals, measured in relatively dust free regions, are very 
similar, implying that stellar populations in bulges and inner disks 
are similar. Temdrup et al. (1994) find the same effect, but with 
more scatter in the colours, and with a sample that includes redder 
galaxies. As we do here, they claim that these galaxies are red be- 
cause of internal extinction. The fact that bulge and disk colours, 
even for these red galaxies, are similar shows that similar amounts 
of extinction are present in both components. If galaxies are seen 
close to face-on, i.e. when most of the light comes from regions 
close to the symmetry plane of the galaxy, the colors of bulge and 
inner disk region both are affected by similar amounts of extinc- 
tion. Images of edge-on galaxies show that the extinction is mainly 
situated in the plane. 

Carollo et al. (2001) present a study of V — H colours in a 
sample of bulges, also from HST imaging. To convert from her 
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Figure 8. Left: Central V-H colour as a function of morphological type. Right: Central extinction Ey_j] as a function of moiphological type. 



V — H AB colours to the Johnson system we use here we need 
to add 1.45 mag. By masking out patches of dust and star forma- 
tion regions, they obtain V — H colours that are blue enough that 
they can be explained by stellar population models, for both types 
of bulges discussed. However, the nuclei have often much redder 
colours, see their Fig. 3, that need considerable extinction. The 
problem with selecting clean regions by hand is that this method 
is subjective and probably not reproducible, giving results that 
are biased towards blue colours. The same is done in Carollo et 
al. (2007), where only one galaxy is found with optical-infrared 
colours redder than those of ellipticals and early-type bulges. 

Dust can have important consequences for the interpretation of 
mass profiles from rotation curves and optical photometry. Without 
correcting for extinction, which is often impossible, the resulting 
stellar M/L ratios will be too large, and show a much larger scatter 
than is due purely by stellar populations. Since the dust extinction 
is largest in the central regions, and HI observations often have a 
limited angular resolution, the effect is most important when in- 
terpreting optical rotation curves. Palunas & Williams (2000) find 
that if they correct the surface brightness by Ai^int = l.Q\og[^) 
the resulting M/L ratios do not vary any more as a function of in- 
clination. In the l/-band the correction should be larger by a factor 
~1.5. Such corrections in general are not very meaningful, since 
extinction varies as a function of radius. If one wants to obtain ac- 
curacte M/L ratios in the centers of galaxies, one either has to work 
in the NIR only, or make detailed extinction corrections. 



standing of the formation of these galaxies. Only the most notewor- 
thy relations are shown here; additional relations, including those 
between disk and bulge parameters, are shown in Appendix B. 

We concentrate, in particular, on how late-type spiral galax- 
ies differ from earlier types. We do this by comparing our galaxies 
with the large literature compilation of Graham & Worley (2008), 
and with the recent sample of Fisher & Drory (2008). Graham 
& Worley present a compilation of about 10 samples, mostly in 
the near-infrared. Fisher & Drory analyse a sample of HST data, 
mostly in the optical. For the relations with central velocity disper- 
sion we compare with the sample of early-type galaxies of Balcells 
et al. (2007b), hereinafter B07. For the literature samples, we took 
the numerical values from the tables available in the published pa- 
pers, adjusting the scale lengths and effective radii in order to ho- 
mogenise the adopted value for the Hubble constant to Hp = 70 
km s ~^ Mpc~^ (for consistency with what done bv iGanda et al.l 
l2006l and with the luminosities calculated in Section HI . The to- 
tal magnitudes (of bulge and disk) of iBalcells et al.l ( l2007bi) were 
converted from K to H using the typical colour for disk galax- 
ies — A" = 0.21 mag (Frogel et al. 1978). For most of the 
relations that we will investigate and show graphically, we list in 
Table [5] the linear Pearson correlation coefficient c, calculated for 
the galaxies of our sample only. For the strongest correlations, with 
c| ^ 0.5, we performed an outliers-resistant linear fit, assuming a 
dependence of the form y — a + bxx,y and x being the variables 
involved and a and b the fit parameters. Those fits, for which the 
parameters are given in Table [6] are discussed in detail below. 



6 BULGE AND DISK PARAMETERS 

Here we study the structural parameters determined in the previous 
Section together with data from the literature. 

We quantitatively investigate the structural parameters ob- 
tained from the bulge-disk decomposition as a function of various 
fundamental galaxy parameters, such as morphological type, cen- 
tral velocity dispersion, and luminosity, to obtain a better under- 



6.1 Trends among the structural parameters 

In Fig. |9]we show a number of structural parameter as a function 
of stellar velocity dispersion. Here one should note that due to the 
paucity of velocity dispersions in the literature we can only com- 
pare 2 samples of about 20 galaxies, of resp. early and late type 
spirals. There are some interesting things to note. Several papers, 
e.g. B07, Fisher & Drory (2008) and Graham & Worley (2008), 
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show that the Sersic index n correlates reasonably well with bulge 
luminosity, with galaxy type and bulge to disk ratio, but with con- 
siderable scatter. Fig. shows that the same holds for the cen- 
tral velocity dispersion a. There are late-type spirals with velocity 
dispersions larger than 100 km/s, that show n-values considerably 
larger than 1. If all bulges of late-type spiral galaxies would be 
disks (also called pseudo-bulges), their n-values would be close 
to 1. This is however not the case, and it appears that the central 
velocity dispersion of the galaxy is determining the shape of the 
surface brightness profile of the bulge, rather than the galaxy type, 
which is determined more by the total star formation or the bulge to 
disk ratio. Consistent with this picture is the relation between bulge 
luminosity and central velocity dispersion. Already B07 found a 
tight relation between these 2 quantities for early-type spiral galax- 
ies. Here we see that late-type spirals follow the same relation. So, 
even though the surface brightness of bulges of late-type spirals is 
lower than of early-type spirals, and the B/D ratio for a given ve- 
locity dispersion is lower as well (Fig. |9j3c), the bulge luminosity 
- (J relation is the same for both samples of galaxies. Since there 
is a strong correlation between the bulge luminosity, or central sur- 
face brightness, or Sersic index n (Graham & Driver 2007), and 
black hole mass, a relation which also holds for the luminosity and 
black hole mass of elliptical galaxies and SO's, it is tempting to as- 
sume that bulges of late-type spirals also lie on this relation, and 
that the formation of those bulges is closely linked to the black 
holes in late-type spirals. If a central black hole is associated with 
every bulge, it would also explain why the tight black hole mass 
vs. central velocity dispersion relation (Ferrarese & Merritt 2000, 
Gebhardt et al. 2000) can be reproduced with bulge luminosity re- 
placing central velocity dispersion, but not with disk or total galaxy 
luminosity. The lack of correlation between the sizes of bulges and 
disk with central velocity dispersion for late type spirals might also 
be a result of the fact that bulges of late-type spirals are less con- 
centrated objects with n ~ 1, which do not follow the homology 
of elliptical galaxies and elliptical-like bulges. 

We also like to point out that it appears that bulges of early 
and late-type galaxies show systematic differences: on the aver- 
age early type spirals with low central velocity dispersion have 
much more concentrated bulges and disks: not only are disks and 
bulges smaller, but also the surface brightness of both components 
is higher. 

6.2 Trends with Stellar Populations 

This sample is unique in the sense that its central stellar popula- 
tions have been studied using stellar absorption lines, which are 
almost independent of extinction, by Ganda et al. (2007), so that 
we can investigate whether the stellar populations depend on the 
structural parameters discussed above. In other words, are the age 
and star formation history of a galaxy related to the way it is built, 
and viceversa? We try to address this issue by plotting the results 
of our bulge-disk decomposition against population parameters. A 
few noteworthy relations are discussed here. The rest of the rela- 
tions are discussed in Appendix B. 

Fig.[TO]shows the Sersic parameter n and the disk central sur- 
face brightness and as a function of the central H/9 absorption, a 
measure of the age of the galaxy, and the star formation time-scale 
r . The stellar population parameters refer to a central aperture of ra- 
dius l'.'5. One sees that n increases of decreasing H/3. This implies 
that for objects with relatively young stellar populations, bulges are 
exponential, and that old bulges are more similar to objects with 
a de Vaucouleurs distribution. This would be the case if the ob- 



jects in which recently stars have formed are all disk-like, i.e. pseu- 
dobulges, and that old bulges are more spheroidal-like. The latter 
statement is very well applicable to bulges of early-type galaxies 
(e.g. Kormendy & Kennicutt 2004, Fisher & Drory 2008), but is 
not so obvious for late-type spirals. It implies that some late-type 
spirals have old bulges, that do not follow exponential distributions. 
Maybe our Milky Way galaxy has a bulge similar to those. We also 
see that n decreases for increasing star formation time scale r. 
Late-type spiral disks, like the disk of our Milky Way, have long 
star formation time scales. As a result, parameters such as H/3, 
which are strongly influenced by the last burst of star formation, 
generally show large values, indicative of the presence of young 
stellar populations. We find that galaxies with small star formation 
time scales, i.e. elliptical-like objects, have high rz-values. Inter- 
esting are also the correlations with disk scale length. If there are 
signs of younger populations in the center, or if the star formation 
time scale is long, then the disks are generally smaller. This effect 
is probably due to the fact that more massive galaxies generally are 
of earlier-type, having a shorter star formation time scale, and show 
older bulges. 

6.3 Trends between structural parameters 

Coming to the relations among the structural parameters them- 
selves, in Fig. [TT] we plot the size and surface brightness of the 
bulges and disks of our galaxies against each other, and add the 
samples of GW and Fisher & Drory (2008) on the right: in the left 
panel we plot the central surface brightness of the disk against the 
effective surface brightness of the bulge; in the right one, disk scale 
length against bulge effective radius. We see a strong trend that 
galaxies with brighter disks harbour brighter bulges, and a tendency 
for galaxies with larger disks to have larger bulges (see before). The 
figure shows conclusively that bulge and disk are related, and that 
e.g. high surface brightness bulges do not co-exist with low surface 
brightness disks and vice-versa. The figure suggests that disks do 
not form without affecting the bulge. If a disk is formed in a galaxy 
which already contains a bulge, its size and surface brightness is ei- 
ther determined by the mass distribution of the bulge, or the bulge 
adapts itself, implying that stars are being formed in this process. 
The strong correlation between ^le, Bulge and iJ.o,Disk seems to sup- 
port the process of the formation of bulges through secular evolu- 
tion of disks. 



7 CENTRAL COMPONENTS 

In recent years, several studies have shown that many spiral galax- 
ies host in their innermost regions an extra component, distinguish- 
able in the images and/or photometric profiles as an excess light 
above the (exponential) disk plus (Sersic) b ulge. Nuclei are com- 
mon both in early- and late- type galaxies. IPhillips et al.l j 1 9961) 
found unresolved bright nuclei in six out of 10 late-type spirals; 
Matthews & Gallagher ( 1997) noticed a compact star cluster in the 
centres of 10 very late-type spirals, within a p rogram studying 49 
objects; out of a sample of 77 Scd-Sm galaxies, iBoker et aH ^20021) 
detected a distinct component in 59 cases. In some of the late- 
type spirals, these clusters are the only prominent source within 
a few kiloparsecs from the centre. In earlier-types their detection 
is c omplicated by the pr esence of a luminous and extended bulge, 
but ICarolIo et al.l ( 11998 ') showed t hat nuclear clusters ar e present 
also there. As already mentioned, iBalcells et al.l i2007ah found a 
central light excess in 90 percent of a sample of SO-Sbc galaxies 
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Figure 9. Structural parameters as a fun ction of central stellar velocity dispersion for our late-type sample (black filled circl es), compared to th e sample of 
early-type spirals of lBalcells et aljE007bl (green asteriscs). The velocity dispersions for our sa mple are taken from T able 1 in lGanda et alj bOOd) . As for the 
morphological T-type, we have artificially put all the galaxies earlier than T = in the sample of lBalcells et alj2007bl to T = 0. 



using NICMOS imaging, detecting both unresolved and extended 
sources; according to the authors, the extended inner components 
are geometrically flat systems and could be inner disks, rings or 
bars; the um'esolved ones are most likely star clusters and are found 
in more than one third of the galaxies with inner component. 

In our sample, two galaxies out of 18 (NGC864 and 
NGC 4102) present a central depression and the remaining 16 show 
a central light excess above the best-fit to the bulge; one of these 
(NGC 5678) is an ambiguous case, possibly because of the very 



extended dust lanes affecting this galaxy, but we included it in the 
'c entral component' gro up. The numbers are very similar to those 
of lBalcells etZl i l2007a '). Out of these 16, eight belong to the sam- 
ple studied by Boker et al. (2002): all the Scd and Sd galaxies in 
our sample. The authors detect nuclear star clusters in all of them 
and publish an estimate of their size and luminosity in the I— band. 

For all the 16 galaxies, we built model images of the bulges 
on the basis of the fit parameters, and subtracted them from the 
HST 'bulge images' (NICMOS or WFPC2-F814W depending on 
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Figure 10. Correlations between the structural parameters n and disk scale length h as a function of stellar population parameters H/3 (in A) and star formation 
timescale r (in Gyr). The dashed lines overplotted represent a hnear fit to our data. 



the object). The images that we obtained contain by construction 
the residuals from the global (exponential) disk + (Sersic) bulge 
fit to the galaxy: dust features, spiral structure, bars, and the inner 
component in which we are interested. We estimated the HWHM 
of the central component using the IRAF rimexam task, resulting in 
values ranging between ~ 0.07 and 0.23". We then calculated the 
flux enclosed in a circle of radius twice the HWHM and converted it 
to a magnitude scale {H~ band). For NGC 4102 and NGC 864 we 
calculated the flux by measuring all the light in the inner aperture 
of diameter 0.227". The computed magnitudes (both apparent and 
absolute) are listed in Table[7] where an asterisk marks the galaxies 



in common with the sample analysed bv lB6keretai]j2002l) . For all 
galaxies for which both HST F606W and F160W data are available, 
we also calculated the central V—H colour in an aperture of diame- 
ter 0.227" (0.500" for NGC 5585), and also tabulated these in Table 
m Fig.[T2]shows the comparison between t he apparent magn itudes 
that we calculated and those published bv lBoker et al.l ( 1200 2) . for 
the inner component of the galaxies in common. The magnitudes 
of Boker et al. refer to the /— band, while our data are H— band 
data (see Table |4l(. The only galaxy in common with Boker et al. 
(2002) is NGC 5585, which is representative for the galaxies in this 
group, and clearly shows much less extinction in itsV — H color 
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map. We also see that the V — H colors of the nuclei are similar 
within 0.2 mag to the colors in the inner aperture with diameter 2.4" 
(A{(V - H)^- {V - H)r^uc) = -0.007 ±0.199 mag). Even 
though only about 10% of the light in the larger aperture comes 
from the nuclear cluster, it shows that the stellar populations of 
the nuclear cluster are not decoupled from those in the rest of the 
galaxy, as is e.g. the case in globular clusters in galaxies. 

Walcher et al. (2006) took high resolution spectra for a number 
of the Scd, Sd and Sm galaxies of Boker et al. (2002) and analysed 
them using stellar population synthesis, allowing the extinction to 
be fitted as a free parameter. They find an average extinction in / 
of ~ 0.4 mag, which is very comparable to what we find (Ay = 
0.68, with a scatter of 0.55): When using the Galactic extinction 
law (Rieke & Lebofsky 1985) the extinction in V is 1.64 x the 
extinction in Ic, or 0.66 mag. Given the fact that the galaxies of 
Walcher et al. are later type than ours, we expect less extinction than 
in our sample. However, large variations probably exist from galaxy 
to galaxy, and also the amount of extinction that one obtains varies 
with method. We can only conclude that extinction is important in 
late-type galaxies. Note that NGC 1042, in common with our paper, 
has by far the largest extinction (A/ = 1.35 mag), consistent with 
Figure [72l 
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8 SUMMARY AND CONCLUSIONS 

In this Paper we performed a photometric bulge-disk decomposi- 
tion for our sample of 18 late-type spiral galaxies. We used NIR 
(H- band) archive images (DSS, 2MASS and HST) from which 
we retrieved photometric profiles covering the whole disk and hav- 
ing a high spatial resolution in the inner parts. We fitted an expo- 
nential disk and a Sersic bulge, using a method that allows the two 



Figure 1 2. Apparent magnit ude of the inner component against those pub- 
lished bv lBoker et alj )2002l) . The solid line overplotted indicates the / = 
H + 1.5 line. The clear outlier below the solid line represents NGC 1042. 

components to have different intrinsic shapes, but without going to 
a full two-dimensional fitting of the light distribution; we studied 
in detail the correlations between the fit parameters, and we inves- 
tigated the high resolution V — H colours in this sample, in the 
nuclear cluster and outside it. From the V — H colour maps and 
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Table 5. Correlation coeiScient c for the parameter-parameter relations, for 
our sample only. We report the correlation coefficient for all the relations of 
Figs. IBH.al and b, and lBE. al and b; for those of Figs. lBBI we only list those 
with |c| 0.5. Columns (1), (2) and (4), (5): variables; columns (3) and 
(6): correlation coefficient. For the relations involving the T-type, all the 
individual points are considered and not the average values within a type 
(as done instead in Fig. IBIl. al and b). 
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log(r) 


n 


2.638 


-1.590 


log 


r) 


log{h) 


0.782 


-0.481 


Mi3„ige 


log((T) 


0.190 


-0.081 


Mb 




T 


15.442 


0.506 


Mc,6 


M0,d 


8.092 


0.527 


log 
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Table 6. Pai'ameters of the Unear best-fit for the relations with |c| ^ 0.5 in 
Table|5] The relation y = a + h y. x, with y and x variables (columns (1), 
(2), (5), (6)) and a and b fit parameters (columns (3), (4), (7), (8)) is fitted 
to the data for our sample only. For the relations involving the T-type, all 
the individual points are considered and not the average values for the same 
type (as done instead in Fig. lBLal and b). 
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3.271 
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15.77 
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Table 7. Column (1): NGC identifier of the galaxies with measurable inner 
component; columns (2); apparent magnitude mi„„er of the inner com- 
ponent, in H— band; columns (3); absolute magnitude yiinner of the in- 
ner component, in H— band. Columns (4); Central colours in apertures of 
0.227" diameter (0.500" for NGC 5585) for the galaxies with both HST 
F606W and F160W data. The asterisks marks the galaxies in common with 
the sample of .Boker et al.. c2Q0Z) . 

the Mg b maps of Ganda et al. (2007) we made extinction maps for 
10 of the galaxies. We summarize below the main results. 

• Late-type spirals can be fitted well by a model consisting of a 
Sersic bulge, an exponential disk, and a nuclear light access. Out- 
side the central regions, the bulk of the galaxy's light is well de- 
scribed by one exponential disk, although in some cases a double 
ex ponential (inner outer d isk) would provide better fits, as tested 
bv lPohlen & Truiillol l liooel) . In most cases bulges are rounder than 
their disks, but for the latest types the opposite sometimes holds, 
possibly in relation with the presence of a bar. 

• Bulges of late-type spiral galaxies follow the same relation be- 
tween central velocity dispersion and bulge luminosity as bulges of 
early-type spirals. This might indicate that the black hole - sigma 
relation also holds for late-type spiral bulges. Late-type bulges 
have, however, lower surface brightness and are larger than early- 
type bulges. 

• The structural parameters of our bulges agree statistically with 
the large compilation of Graham & Worley (2008), showing a 
strong correlation between surface brightness of bulge and disk, 
and a less strong relation between bulge effective radius and disk 
scale length. 

• The star formation time scale inversely correlates with the 
Sersic index n, implying that exponential bulges have long star 
formation time scales, in agreement with the prediction from sec- 
ular evolution models, that exponential bulges are disks forming 
through secular evolution. 

• In 16 out of 18 galaxies the galaxy's profile presents an excess 
light with respect to the fitted disk -I- bulge, which can be inter- 
preted as an additional tiny component, in many cases (in all the 
later types) a nuclear star cluster. We give an estimate of the mag- 
nitude of these inner components. 

• Using V — H colour maps from HST imaging and SAURON 
Mg b maps, we have determined model-independent extinction 
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maps for 10 of the 18 galaxies. We find that the central Ay ranges 
from to 2 mag, with many galaxies being optically thick in the 
optical. The V — H profiles show a lot of structure, mostly due to 
extinction, with much larger gradients than are generally displayed 
by early-type galaxies. 

• The colours of nuclear clusters are the same as the inner re- 
gions of their host galaxies, with a scatter of 0.2 mag, indicating a 
similar composition and age as the stellar populations just outside 
the center, unlike e.g. galactic globular clusters. 
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APPENDIX A: INTERNAL CONSISTENCY OF THE 
METHOD - LITERATURE COMPARISON 

In order to ensure the reliability of the method applied, we searched 
the literature for published H— band profiles. For the late-type 
spiral galaxies of our sample extremely little is available. One 
of the few existing references is the work of IJunt & M alkan 
( l2004h . who present NICMOS photometry for 250 active and nor- 
mal galaxies, among which six are in common with our sam- 
ple: NGC 488, 772, 2964, 3949, 4102, 5678. In almost all cases they 
used the same images as we did. The authors were so kind to pro- 
vide us in tabular format the profiles for those galaxies, that they 
extracted keeping the centre of the isophotes fixed, but allowing 
the geometric parameters to vary over the whole radial range. In 
order to match their approach, we extracted again photometric pro- 
files from the 2MASS and NICMOS images (DSS profiles are here 
redundant, since the literature profiles are available only in the in- 
ner ~ 10" only; therefore we did not use them) a long the average 
positio n angle and ellipticity of the profiles from lHunt & MaUcarj 
( l2004h. InFi g. lAT] we show for those six galaxies the comparison 
with lHunt & MalkanI <2004l) : the profiles (left panels) and the dif- 
ference between them, plotted on the common radial range. The 
agreement is satisfactory (see also figure caption) . NICMOS pro- 
files f or five of the six galaxies in common with ' Hunt & MalkMil 
( l2004h (NGC 488, 772, 2964, 3949, 5678), again deri ved from the 
same images that we used, are published also by ISeigar et al.l 
( l2002h : from a visual inspection of their profiles, we noticed that 
there is a significant offset with ours and Hunt's profiles, which 
is probably due to the adoption of another magnitude system (AB 



magnitude system), therefore we did not proceed to a more detailed 
comparison. 

As a reference for the groun d-based H— b and photometry, we 
chose to use the work of Gava zzi et al.l ( 1200 ll) . who present H— 
band observations and surface brightness profile decomposition for 
75 faint galaxies, mostly dwarf ellipticals belonging to the Virgo 
cluster. They collected H— band images at the ESO-NTT telescope 
in La Silla (Chile) and at the TNG on La Palma (Spain), reduced 
and calibrated the images and extracted brightness profiles fitting 
ellipses with ellipse centres, ellipticities and position angles left as 
free parameters. This corresponds to the first step in our procedure, 
as we described above. Since there is no overlap between their sam- 
ple and ours, we chose some galaxies out of their paper (NGC 4353, 
4440, 4600, 4706 and 4743), retrieved the 2MASS and DSS images 
for these galaxies and extracted brightness profiles with free geo- 
metric parameters, subtracted an estimate for the sky background, 
combined them into a single profile and converted the latter to an 
absolute magnitude scale in the same way as done for our own 
galaxies -we only skipped the steps where we fix the geometric pa- 
ramete rs, for consistency with the approach used bv lCavazzi et al.l 
( l200l[) . Fig. IA2I shows for the chosen galaxi es the compar i son be - 
tween our profiles and those published by iGavazzi et al.l ( 1200 ll) . 
kindly made available in tabular format by the authors. The figure 
shows that the agreement is good, confirming the reliability of our 
method. 



APPENDIX B: STRUCTURAL PARAMETER RELATIONS 

In this appendix we present a large number of relations, between 
structural parameters and luminosity of bulge, disk and total galaxy, 
and between structural and stellar population parameters. A few of 
these have already been discussed in Section 5. For the relations be- 
tween structural parameters and luminosity of bulge, disk and total 
galaxy we find that our late type spiral galaxies generally behave in 
the same way as the compilation of Graham Sl Worley (2008). In 
Section 5.1 and 5.2 the most important correlations were discussed. 
Here we discuss more relations, for completeness. 

We start looking at the relations between the structural param- 
eters of bulge and disk and global galaxy properties such as mor- 
phological type and velocity dispersion. We plot these parameters 
against each other in Fig. IB 1. al and IB 11 .51 for our sample with the 
compilation of GW always plotted as well to guide the eye. If we 
focus on our sample only, we can notice some trends: n tends to 
decrease going to later types, and both the effective surface bright- 
ness of the bulge and the central surface brightness of the disk de- 
crease going to later types. These correlations agree with GW and 
B07. Interesting to note is that for galaxies earlier than Sbc sur- 
face brightness of bulge and disk is virtually independent of type. 
For later types there is a strong correlation. There galaxies become 
smaller, and have lower surface brightness of both bulge and disk, 
for increasing morphological type (see also de Jong 1996, paper 
III). We confirm that there is no strong corelation between the size 
of bulges and disk and morphological type. Also, the size ratio r^lh 
is independent of type, and we do not find the mild decrease that 
iMacArthur et alj ( l2003h claim. See GW for an extended discussion 
about this ratio. The most important correlations with central celoc- 
ity dispersion have been discussed in Section 5.1. We find that the 
effective radius of the bulge does not correlate with a, nor does the 
scale length of the disk, or the ratio of bulge effective radius to disk 
scale length. There is a good correlation between bulge luminosity 
and a (see also B07) and between bulge to disk ratio and a. 
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Figure Al. Comparison between our combined NICMOS + 2MASS profiles and the profiles published by iHunt & MalkanI j2004l ) for the six galaxies in 
common, plotted on the common radial range. In the left panels, the crosses represent our profiles, and the filled triangles those kindly made available by 
Leslie Hunt, also plotted on the common radial range. In the right panels, we plot the difference between the two. The difference is relevant only in the 
innermost 0.1-0.3". and might be due to small differences in the centering and/or in the geometric parameters. In particular, NGC3949 represents the case 
with the largest difference in the innermost points, and it is also the case where the difference between the elhpticity profile from Hunt and the mean value we 
adopted is the largest. 



In Fig. IBE.al and IBE.bl we plot the structural parameters of 
bulges and disks as a function of bulge, disk and total Jf-band lumi- 
nosity. We see that surface brightness of bulges and disks increases 
with increasing bulge or disk luminosity, consequence of the con- 
stancy of disk scale length and bulge effective radius. Since later 
type spiral galaxies tend to be fainter, we find that the shape pa- 
rameter n increases with the luminosity of the bulge, disk or the 
whole galaxy. The relation between the bulge magnitude Mb and 
n (top right panel in Fig. lBE.aT i is particularly tight: the correlation 
coefficient c is ~ —0.7, and it becomes even tighter (c ~ —0.92) 
when removing the outlier point representing NGC 864. As shown 
by B07 in their Figure 6c, spiral galaxies follow the same trend 
as the Virgo ellipticals from lCaon. Capaccioli & D'Onofriol ( Il993h 
in the plane {Msuige, n). We find a marginally significant corre- 
lation between bulge luminosity and effective radius, in the sense 
that brighter bulges are also larger. When the whole compilation of 
GW is considered this correlation is not found. In general, when 
the bulge luminosity increases, the Sersic index n decreases, and 
the surface brightness increases, with constant effective radius. For 



late type spirals only, which all have a similar Sersic index (n=l), 
there is apparently a tendency that both fie and re increase. The 
marginally significant correlation between disk luminosity and disk 
scale length is probably caused by the two bright outliers. The B/D 
luminosity ratio increases as galaxy bulges become more lumi- 
nous, while it does not show strong trends with the disk magnitude 
Moisk- In general, from this work we infer that several correlations 
between structural properties of spiral galaxies are tighter versus 
M Bulge than versus Mntsk (see T able [5]l, in agreement with the 
conclusion o f iBalcellsetai] ( 12007 1/). This is particularly intriguing, 
implying that (some of) the properties are driven by the mass of the 
bulge. GW, who find the same strong relation, also comment that 
these new decompositions, with Sersic rather than de Vaucouleurs 
bulge, show that very few spirals have B/T > 1/3. In a long dis- 
cussion GW confirm this conclusion, and point out that the B/D 
ratios in cosmological galaxy simulations are on the average much 
larger, showing that the models, or probably the recipes for super- 
nova feedback, need to be improved considerably. 

Fig. lBBI shows the Sersic parameter n, the bulge effective sur- 
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Figure A2. Comparison between the combined 2MASS-DSS profiles and the profiles from lGavazzi et al" for five galaxies taken from Gavazzi et al. 

(2001, see text for a more detailed description). In the left panels, the crosses represent our profiles, the filled triangles those from Gavazzi et al. (2001); the 
vertical scale is i n absolute H magn itudes arcsec"^ . In the right panels, we show on the common radial range the difference in magnitudes arcsec"'^ between 
the profiles from lGavazzi et al]i2001.) and ours. In most cases, the differences are significant (more than 0.1 mag arcsec"^) only in the innermost 2-3", where 
our 2MASS profiles are affected by seeing smearing m uch more than the ones of lGavazzi et alj j200lh . and at radii larger than 25-30", coiresponding to the 
external limit of the profiles from lGavazzi et aU i200lh . This proves that the agreement is good. 



face brightness and radius, the disk central surface brightness and 
scale length, the r^/h ratio and the B/D luminosity ratio as a func- 
tion of H/3, AMgfc', age and star formation time-scale r. The stel- 
lar population parameters refer to a central aperture of radius l'.'5. 
AMgfc' is the difference in magnitudes between the observed Mgb 
index and the prediction, based on the obser ved velocity dis persion , 
given by the Mg2 — <j relation published bv lj0rgensen et alJ jl996h . 
which is now believ ed to hold for old spheroids (see discussion in 
|Peletieretal1 ( l2007h '): AMgi' is therefore a measure of the impor- 
tance of young populations in the galaxy. For the sake of clarity, 
we remind the reader that AMgZ?' can be derived from Mg b and 
a in the following way, adopting the Mg2 — a relation from Jor - 
gensen, expressed as in Eq. 3 in the paper of lGandaetai] ( l2007h . 
the definition of Mgb' inferred from Eq. 2 in the same paper and 
the conversion between Mg2 and Mgb' given in Eq. 4 there: 

AMgfc' = -2.5xlog(l-Mgfc/32.5)-0.096xlog(cr)+0.062.(Bl) 

Going back to Fig. IBB! and to the description of the quanti- 
ties involved here, the age is the Single Stellar Population (SSP)- 
equivalent age inferred by comparison of observed line-strength in- 



dices with models and r is the e-folding time-scale for star forma- 
tion, assuming a galaxy age of 10 Gyr and an exponentially declin- 
ing star formation rate with tim e. For details we re fer the reader 
to Sections 5.3, 6.1 and 6.3, in iGanda et"al] ( |2007|) . We can ob- 
serve some trends, a few of which significant (|c| 0.5), in those 
cases we have also plotted the correlation as a dashed line: when 
H/9 increases, i.e. when the stellar population becomes younger, n 
decreases, both /icDisfe and fie, Bulge become fainter, and the disk 
scale length h and the bulge effective radius Ve decrease. There are 
also indications that both Ve and h become smaller with AMgfe' 
becoming more negative, which means that galaxies farther away 
from the Mg2 — o relation of old ellipticals, i.e. younger galaxies, 
have both smaller bulges and smaller disks. No relevant trend is 
found with age: possibly a confirmation of the fact that, in general, 
SSP-equivalent ages and metallicities a re not a good approx ima- 
tion for spir a l gala xies, as discussed by jPeletier et al.l ( |2007|) and 
iGanda et al.l ( l2007h . A more relevant parameter for our galaxies 
seems to be r, in the framework of a continuous star formation: 
as r increases (and the star formation history is less and less well 
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described by a single burst), n becomes smaller, /io_d and /^e^;, be- 
come fainter and and h become smaller. No significant trend is 
observed involving the ratio re/h ratio and the B/D luminosity ra- 
tio. Globally, Fig. lBBI seem to corroborate a scenario where galax- 
ies with a quiescent, almost constant star formation, hosting com- 
posite stellar populations, with (at least) some young stars, have 
almost-exponential bulges and small and faint bulges and disks, 
while galaxies that formed stars in a single burst long ago (more 
SSP-like) tend to have higher n and bigger and brighter bulges and 
disks. A similar scenario, in which the character of the star for- 
mation history is rep laced by the galaxy age, is suggested also by 
iHunt. Pierini & Giov anardi (2004), as we mentioned earlier: expo- 
nential bulges would be young products of secular evolution; as the 
galaxy ages, the bulges would grow and their n increase. 
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Figure O. a. Bulge and disk structural parameters as a function of morphological T-type and central stellar velocity dispersion cr. The black filled circles 
represent our sample; blue tiiangles represent the sample of Fisher & Drory (2008). The red open squares represent the literature compilation of Graham & 
Worley (2008). Left panels: from top to bottom we plot n, fie, Bulges ''e and fJ.o,Disk'^ for s^ch sample and each value of T-type, we plot the average value 
of the quantity of interest for all the galaxies of that type in that sample, with an errorbar indicating the standard deviation aro und the mean . Right panels: we 
plot the same quantities against log{a), measured in km s^ ^ . The velocity dispersions for our sample are taken fr om T able 1 in lGanda et alj i2006l) . As for the 
morphological T-type, we have artificially put all the galaxies earlier than T = in the sample of .Balcells et al.2007bl to T = 0. 



The nature of late-type spiral galaxies: structural parameters, colour profiles, and dust extinction 25 



I* IB 



15 1.5 3 



1 y Li LO 



i.'j i.ii d y.y 



4 S 6 LO 



IS l.E- 3 



I ' ' ' I ' ' ' I 



■1 U Li LO 

T- Lyp^- 



i.'ii L.Li ^ y.y 

at; .j(Ltm/n| 



Figure|B]l.b. Bulge and disk structural parameters as a function of morphological T-type and central stellar velocity dispersion a (see Fig. lBH.a) . Left panels: 
from top to bottom we plot disk scale length, bulge effective radius to disk scale length ratio, //-band bulge-to-disk ratio and bulge luminosity against T; right 
panels: same quantities plotted against \og(u)m km s^^. 
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Figure|Bf .a. Correlations between structural parameters of bulge and disk with the luminosity of the disk (left), bulge (middle) and the whole galaxy (right). 
For symbols see Fig. IBt.a1 From top to bottom we plot coirelations of Sersic index n, bulge effective surface brightness, bulge effective radius and central 
disk surface brightness. The sample of Graham & Worley has been plotted in different colours, ranging from red for the earliest type galaxies to blue for the 
latest types. 
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Figure|Bf .b. Correlations between structural parameters of bulge and disk with the luminosity of the disk (left), bulge (middle) and the whole galaxy (right). 
Left panels: from top to bottom we plot correlations with disk scale length, ratio of bulge effective radius to disk scale length, and //-band bulge-to-disk ratio. 
The sample of Graham & Worley has been plotted in different colours, ranging from red for the earliest type galaxies to blue for the latest types. 
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Figure|Bf. Correlations between the structural parameters n, /.te ^, re, /^o Disfc. h, re/h and the B/D luminosity ratio against the population parameters H/3 
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